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We prove basic theorems about the ground states of the S — 1 Bose-Hubbard model. The results 
are quite universal and depend only on the coefficient U2 of the spin-dependent interaction. We 
show that the ground state exhibits saturated ferromagnetism if U2 < 0, is spin-singlet if U2 > 0, 
and exhibits "SU(3)-ferromagnetism" if U2 = 0, and completely determine the degeneracy in each 
region. 
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Recent progress in cold atom experiments, in partic- 
ular those of bosons trapped in an optical lattice, has 
led to a renewed interest in the low energy properties of 
the Bose-Hubbard model [lJ-OI . Especially the system of 
spinor bosons, in which hyperfine spin degrees of freedom 
couple to many-body quantum physics of bosons, are ex- 
pected to have a variety of phases including the spin- 
singlet, nematic, and ferromagnetic phases [4J48J. From 
a theoretical point of view, it is also interesting to com- 
pare the situation with that of the Fermi-Hubbard model, 
where intricate interplay between the spin degrees of free- 
dom and many-body quantum physics has been studied 
in depth 

Here we shall present basic theorems about the ground 
states of the standard Bose-Hubbard model for 5 = 1 
bosons. We precisely characterize, in each parameter 
range, the degeneracy and the total spin angular mo- 
mentum of the ground states, and identify the signs of 
the "wave function" represented in an appropriate basis. 
The result depends only on the single parameter U2, and 
one gets ferromagnetic ground states for JJ% < 0, a spin- 
singlet ground state for U2 > (when the boson number 
is even), and ground states with "SU(3)-fcrromagnetism" 
for U2 = 0. These results are surprisingly universal, 
and do not depend on the lattice structure, the spin- 
independent interaction, or the boson number. This is in 
a sharp contrast with corresponding results in the Fermi- 
Hubbard model [§-11 1 and the Heisenberg spin system 
12l Il3| , where strict restrictions on lattice structures and 
the electron numbers are usually required. 

We hope that these basic and rigorous results will be 
useful for the future theoretical and numerical studies of 
the Bose-Hubbard model. We note that our Theorem 3 
for U2 = is a special case of the theorem by Eisenberg 
and Lieb proved for spinor bosons in continuum [14) . 

Definitions and main theorems: We consider a sys- 
tem of N spinor bosons with S = 1 on a finite lattice A, 
where N is arbitrary and fixed. For each site x G A, 
we denote by a\ a and a XM the creation and the an- 
nihilation operators, respectively, of a boson at x with 
spin a E {+,—,0}. We define the number operators 
by n Xta := at a x , a and n x := 



the spin operators by S x a ^ ■= J2a.Te{+.-.o} a x,a 



0} n x,cn and 
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for a — x, y, z, where Sa 
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5 = 1. We write S x := {S x * ] , S [ x y \ S ( x ] 
?(<*) 



are the spin matrices for 
The total 

spin operators are S^t '■— ExeA^^ '• ^ s USU& 1 we 
define 5j^ t = S^ot ^ ^toti an d write the eigenvalue of 
(Stot) 2 = E Q=x , y , z (^t ( o Q t ) ) 2 as S tot (S tot + 1). 
We shall make a frequent use of the basis states 
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where $ V ac is the state with no bosons in the trap, and 
the multi-index n = (n x ,a-)xeA, o-ef+,-,0} is a collection 
of nonnegative integers such that n Ij(T = AT. The 

whole Hilbert space H is spanned by all such <£>„. For 
M = 0, ±1, . . . , ±AT, we denote by Im the set of n satis- 



fying Ex( 



-) = M, and by Hm the subspace of 



% spanned by with n £ Im- Note that any $e?i 
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satisfies = M*, and that U = 

We consider the standard Hamiltonian of the 5=1 
Bose-Hubbard model 0, @ 

H = - t x ^ y a x a a v , a + Vx n x 

x,yeA, o-e{+,-,o} 
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which has a global SU(2) symmetry of the spin rotation. 
We assume that the hopping matrix elements t x . y = t y . x 
is nonnegative, and the whole lattice A is connected via 
nonvanishing t x . y . The standard choice, where one sets 
t x , y = t > for neighboring x, y, and t x<y = otherwise, 
is sufficient. Here V x £ K is the on-site single-particle 
potential, and Uq and U2 are real coefficients for the spin- 
independent and spin-dependent two-body interactions, 
respectively. We make no assumptions on V x and Ug. 
Assumptions on U2 are stated in the theorems. 

In an (idealized) experimental realization of the Bose- 
Hubbard model, one initially prepares specified numbers 
of bosons with +, 0, and — spins in the trap, and lets 
the system evolve in time. As long as U2 7^ 0, spin flips 
take place and the state evolves into a superposition of 
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various spin states, but the eigenvalue of S^ t , which is 
M, is conserved. Therefore the system is expected to 
evolve into a low energy state within each subspace % m ■ 
It thus makes sense to consider local ground states within 
each Hm- 

Theorem 1 — If J7 2 < 0, the local ground state ^/ in 
Hm is unique, and is written as 



* — 'nelM 



(3) 



with a n > 0, and has the maximum possible total spin 
Stot = N. The local ground state energy Efp is indepen- 
dent of M. Thus each 9^ is the global ground state in 
W, and the ground states are (2N + l)-fold degenerate. 

In short the ground states exhibit saturated ferromag- 
netism, and are unique apart from the trivial (2 Stot + 1)- 
fold degeneracy. Note that the ground states of the 
present bosonic system with U2 < may be in var- 
ious phases, including Bose-Einstein condensation and 
Mott insulator. As for the magnetic property, however, 
the ground states universally exhibit saturated ferromag- 
netism. 

Theorem 2 — If U2 > 0, the local ground state ^ff in 
Hm is unique, and is written as 



* — 'n£T M 



(4) 



with fi n > 0. The state *f f s has total spin S to t = \M\ if 
TV — M is even, and S tot = \M\ + 1 if N - M is odd. The 
local ground state energy satisfies E M S = E2% f°r an Y 
M, < £GS UN — M is even, and = E<jf |+1 
if N - M is odd. 

Thus is the the unique global ground state in H 
if N is even, while ^f s , ^q S , and $>®f form a triplet 
of global ground states if N is odd. In short the ground 
states have the smallest total spin as possible, as is ex- 
pected from the antiferromagnetic nature of the inter- 
action. We remark, however, that the present theorem 
is not sufficient to determine the magnetic property of 
the ground states, since the constraint on S to t alone does 
not specify magnetic structure uniquely. Let us recall 
that the Heisenberg model on a bipartite lattice satis- 
fies a theorem very similar to Theorem 2 12, 13j], but it 
may exhibit various magnetic properties (e.g., antiferro- 
magnetic long range order, quantum criticality, Haldane 
phase, dimerization) depending on details of the model. 

Finally we focus on the special case U2 = 0, where 
the spin-dependent interaction is absent. In this case no 
spin-flip takes place during the time evolution, and the 
numbers of +, 0, and — spin bosons are conserved sepa- 
rately. We thus have to talk about local ground states in 
smaller subspaces. For three nonnegative integers U, D, 
and Z such that U + Z + D = N, let us denote by 1u,d,z 
the set of n such that J2 X n x,+ = U, J2x n x- = D, and 
J2 X n x .o = Z. We also denote by Hu,d,z the subspace 
spanned by with n g lu d z- 



Theorem 3 — If J7 2 = 0, the local ground state ^% S D z 
in Hu.d.z is unique, and is written as 



GS 

U.D.Z 
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with 7„ > 0. The local ground state energy E^ S D z is 
independent of (U, D, Z). Thus each ^if D z is the global 
ground state in H, and the ground states are (N+1)(N+ 
2) /2-fold degenerate. 

This theorem is a special lattice version of that proved 
in where details about the degeneracy is omitted. 
See also [IH . The high degeneracy of the ground states 
can be understood as a manifestation of the SU(3) sym- 
metry of the model with U2 — 0. The ground state does 
not exhibit magnetic orderings in the standard sense, but 
there certainly is an "exchange interaction" which real- 
izes the property ([5]) . One may say that the ground states 
exhibit "SU(3)-fcrromagnetism." 

Limiting cases: While the above theorems provide 
only general information about the ground states, much 
more detailed information may be obtained in some spe- 
cial cases. 

As a simple example, we show that one can precisely 
characterize the ground state in the "singlet phase" by 
using a rigorous perturbation theory. Let A be a uniform 
lattice, such as the cubic lattice with periodic boundary 
conditions, and denote the number of sites as |A|. We 
consider the case when v := N/\A\ is an even integer. 
We set V x = for all x e A, and t x , y — t > when x 
and y are neighboring sites and t x y = otherwise. Then 
we shall perturb around the trivial model with Uq > 0, 
U2 > 0, and t = 0. By applying the general results (see 
Theorems 4.1 and 4.3) in [16J, we can show the following. 

Theorem 4 — There are positive constants B and 8 
which depend on the nature of A but not on the size 
|A|. When one has Uo/t > B and f/a/i > B, the unique 
ground state is a small perturbation of the product state 
where each site is occupied by v bosons which form a 
spin-spin singlet. There is an energy gap larger than S 
above the ground state energy. 

Let us also make a remark on the weak coupling limit, 
which essentially is a rewording of the earlier results in 
17, Let us start from a model that satisfies the 

conditions for Theorems 1 or 2, and then take the limit 
Uq.,U2 — > 0. Since the limiting theory is certainly that of 
free bosons, the (local) ground state in "Ho is written as 



Vn 



n!(4=2) 
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where we have assumed that N is even. Here n is summed 
over even integers from to A, and b\ is the creation 
operator of the unique single-particle ground state with 
spin a. Interestingly, the fact that $ has Stot — or N 
alone determines the coefficients r\ n uniquely. 
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We have (apart from normalization) r\ n = 
2 n ? 2 /{Vri.{(N _n)/2}\} if Sua = N, and 
r) n = (-2)-"/ 2 \/^!/(n/2)! if S tot =0 0. Note that 
\r] n \ 2 is proportional to the probability of observing n 
spin-0 bosons in the state <J>. If St. t = N, this probability 
distribution is \r] n \ 2 ~ const. exp[-27V{(n/iV) - (1/2)} 2 ] 
for large N. If S to t = 0, on the other hand, it becomes 
the power law distribution |?7 n | 2 ~ const. / \/n, which is 
quite characteristic. 

Of course exactly the same claim applies to a contin- 
uous model of bosons, i.e., bosons in a standard optical 
trap (see the final part of the paper). To our knowl- 
edge the power law distribution of the number of spin-0 
bosons in the ground state with Stot = has not yet been 
observed experimentally. 

Proof of Theorem 1: We shall show below that 
($„., H$ n i) < for any n, n', and that the whole set 
Im are connected via nonvanishing , H<& n > ) . Then 
the Perron-Frobenius theorem [2(| guarantees that the 
ground state in Hm is unique and is written as ([3]). 

To determine the total spin and the degeneracy, note 
first that (which consists only of spin + bosons) 

has S^t = N. With (J3J in mind, one finds for M = 
-N,...,N - 1 that * M := (S^ ot ) N - M ^ s is nonvan- 
ishing and admits the expansion as in © with nonneg- 
ative coefficients. Noting that j3 n > in ([3]), we find 
(^M,^m) 7^ 0- But we know that ^m is an eigenstate 
of H (since [H, S t ~ t ] = 0), and the uniqueness of the 
ground state within T-Lm shows that \I/ m = (const. )^m- 
Obviously ^ M has Stat = N since [(S to t) 2 , S'totl = °- 

It remains to prove the claim about the matrix el- 
ements. Note that the off-diagonal matrix elements 
($„., H$ n i) in a bosonic system can be directly read 
off from the representation of H in terms of and a's. 
This is in a marked contrast with fermionic systems, in 
which intricate behavior of fermionic sign plays an es- 
sential and nontrivial role. Clearly the hopping term 
yields a matrix element —t x _ y which is non-positive. To 
examine the spin-dependent interaction, we recall that 
(Sz) 2 = {S+S- + S~S+)/2 + (S (z) ) 2 , in which the sec- 
ond term only gives diagonal matrix elements of H in the 
<&„ basis. Noting that 5+ = \f2{a x+ a x ,o + a x [) a x -), 

S~ = \/2{a x Q a x ,+ + a x _a x ,o), an inspection shows that 
the only contribution to off-diagonal matrix elements 
comes from C/a{<4 +«4 _(a x ,o) 2 + h.c.}. This has the de- 
sired sign if Z7 2 < 0. 

To show the connectivity, fix a site x. From any 
configuration n, one first brings all the bosons to x. 
Then by applying the U2 term, one can change the spin- 
components to any desired one (within the fixed M). 
This shows that any configuration in Im is connected 
(via nonvanishing matrix elements of H) to the one con- 
figuration where all the particles sit on x. 

Proof of Theorem 2: It is apparent that the U2 terms 
now produce positive matrix elements ($„,,_/?<!>„/), and 



the Perron-Frobenius theorem does not apply as in the 
above proof. Instead we define new basis states by 
<f>' n ■= (-l)£*"*.°/ 2 $ n , and show that (<f>' n ,H<f>' n ,) < 
for any n,n'. To see this first note that the hopping 
term again yields the matrix element —t Xj y which has 
the desired sign, since the hopping preserves J2 X n x,a- 
Next note that the off-diagonal contributions from the 
U2 terms change the number of spin bosons by two. 
Thus the prefactor (— ™»,o/ 2 yields an extra —1, and 
the matrix element becomes —U2, which has the desired 
sign. Since it is obvious (from the proof of Theorem 
1) that the whole set Im are connected via nonvanishing 
(Q'n, H$' n ,), we can apply the Perron-Frobenius theorem 
to conclude that the ground state in Hm is unique and is 
written as ^ . We note that this ar gum ent is the same as 
that used for S — 1 spin systems in [2ll423| , but seems to 



be more straightforward in the present context of many 
boson system. 

Since H, S^l, and (S to t) 2 are simultaneously diago- 
nalizable, the uniqueness of the ground state (in Hm) 
implies that fy^F is an eigenstate of (Stot) 2 - To deter- 
mine the eigenvalue, consider a toy model on the same 
lattice A, with the same boson number N, and the Hamil- 
tonian obtained by setting t XtV = for all x, y G A, and 
V x = for all x G A\{o} in H of where o is a fixed 
site in A. We still have U2 > 0. When V is negative and 
large enough, the toy model has ground states in which 
all the N particles occupy o and are coupled to minimize 
(S ) 2 . The bosonic symmetry of the state implies that 
the minimum possible Stot in Hm is \M\ or \M\ + 1 when 
N — M is even or odd, respectively. On the other hand, 
the Perron-Frobenius theorem can be applied also to the 
toy model to show that the ground state ^ff in %m is 
unique and admits an expansion like Q with j3 n > 0. 
This means (*^ S , *^ S ) 7^ 0, which implies that has 



the same Stot as \& 



GS 
M 



°tot w M+l fc 



u 



By using S t t t S t ~ ot 



tot/ 



1*1 



The symmetry Efp = E_m is obvious. To prove 
the ordering of E M S , take M > and define * := 

tot^tot = (Stot) 2 ^ 
vjjGS 0+ 0- ^gs v _ 

*Af+H tot°tot*M+l/ — 

2(2M + 3)||*^ b +1 || 2 if N - M is even, and ||*|| 2 = 
2(M + 1)||*m+iII 2 XN-M is odd. We thus see * ^ 0. 
Note that * e "Hm is an eigenstate of H with the 
eigenvalue Efp +1 because [S^ ot ,H] = 0. We thus find 
When N - M is even, * (like *m+i) 



pGS 



has 5 to t = M + 2 while has S to t = M. Thus 

1 $>2i' and * are orthogonal, and the uniqueness of the 
local ground state in TLm implies E M S < E M S +1 . When 
N - M is odd, define f := St^ GS 



S^tS+t = (Stot) 2 - (^tot) 2 - Stot and recalling that 



S^af G Hm+i. By using 



/v& us q+ vi/ GS \ 

\*Af ' J tot J tot W M , 



has Stot = M + l, we find ||*'|| : 

2(M + 1)||*^ S || 2 ^ 0. Since G U M +i is an eigenstate 



we find EfP +1 < E M 



GS 



of H with the eigenvalue Eff , 
This, with the previous inequality, implies the desired 
equality Efp = E M S +1 . 
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Proof of Theorem 3: The proof is essentially the same 
as that of Theorem 1. We use the fact that H$> n >) < 
for any n,n' , and that the set Iu,d,z are connected via 
nonvanishing H& n >). Then the Perron- Frobenius 
theorem proves the uniqueness of the local ground state 
in each Hu,D,z an d the property ([5]). To show the de- 
generacy, we again start from the fully polarized ground 



state *gs 0)0 and observe that (S+^-) D (S+^ u ) z ^ N<0 ll 
is a constant times ^% S -n-z D Z- 



duced two lowering operators S. 



tot 



+- s-0\Z, T( GS 
tot ) *JV,0 

Here we have intro- 
: = Ex4,- a x,+ and 



Ez a x o a x,+> wn i cn commute with H if {/ 2 = 0. 



Some extensions: One can replace the interaction 
part in (0) by E^aI^ W™* - l)/2 + u£°{(S x ) 2 /2 - 
n x }], thus making the interaction site dependent. The 

(2) 

conditions for the theorems then read "if u x < for any 

(2) (2) 

x and u x < for some x" for Theorem 1, "if u x > 

(2) 

for any a; and Mi > for some x for Theorem 2, and 



"if u x = for any x" for Theorem 3. ui"' is arbitrary. 

One can take into account non-local and many-body 
interactions to our theory without any modifications. 
The spin-independent interaction Uon x (n x — l)/2 can be 
replaced by any function U ((n x ) X £\) of the number oper- 
ators, and still Theorems 1, 2, and 3 are valid as they are. 
If one replaces the spin-dependent interaction [^(S^) 2 by 
Ex y Jx,y^>x ■ Sj,, Theorem 1 is still valid under the condi- 
tion J XiV < (and J x<y < for some x,y), but Theorem 
2 no longer holds under any conditions. This is because 
S x ■ S y with i / j generates an off-diagonal term like 
a x o a y,+ a x,+ a y,Oi whose sign is not changed by the intro- 



(0) 



duction of the prefactor (— 1)^* 



>,o/2 



When an external magnetic field in the z-direction is 
applied to the system, one should add to H the new 
terms J2 x [~PA n x,+ ~ n x -} + q x {n x . + + u x> _}], which 
are the linear and quadratic Zeeman terms, respectively 
24]. With this modification, which breaks the global 
SU(2) symmetry, we can only prove the basic properties 
([3]), dU, and (O, as well as the uniqueness of the ground 
state in each subspace. 

For the system of bosons with spin S > 2 and the same 
Hamiltonian @, one can prove theorems corresponding 
to Theorems 1 and 3 in the same manner. Now the model 
at U2 = has a higher SU(25 + 1) symmetry, and the de- 
generacy becomes larger. There seems to be no straight- 
forward extension of Theorem 2, s ugg esting richer phase 
structures in higher S models 11, 2M, HJ . 



Continuous systems: Since our theorems cover mod- 
els with a wide range of parameter values, it is expected 
that they are valid in the formal continuum limit. More 
precisely consider a system of N spin-1 bosons in R 3 with 



the standard Hamiltonian [271. |28| 

H = J d 3 rY,[-^Ur)^ a (r) + V(r)^l(r)^ a (r) 

(T 

d 3 r d 3 r' Mr, r')*t(r)*^ (r')* CT (r)^(r') 



(7) 

where the trap potential V(r) and the spin-independent 
interaction Uo(r,r') are completely arbitrary. Then we 
state (i) when U2{r,r') < for any r and U2(r,r') < 
for some r, all the statements of Theorem 1 are valid as 
they are; (ii) when U2(r,r') = c 2 S(r — r') with c 2 > 0, 
all the statements of Theorem 2 are valid as they are. 

One can easily construct "proofs" based on the stan- 
dard variational argument (see, e.g., [3 [2^]) of these 
statements. We shall not, however, claim that we have 
proved them rigorously since mathematically rigorous 
treatment of continuous systems requires some nontriv- 
ial preparations. In fact the claim (i) may be justified 
by standard techniques (see, e.g., [30|) with suitable as- 
sumptions on the potentials. But a rigorous justification 
of the claim (ii) may require nontrivial effort because of 
the delta function interaction. 

Discussions: We have proved basic theorems about 
the ground states of the standard S — 1 Bose-Hubbard 
model. The result is surprisingly sharp, and we were 
able to completely characterize the finite- volume ground 
states for U2 < 0, U2 > 0, and U2 — 0. A key to un- 
derstand this sharpness may be the large degeneracy at 
U2 = 0, which reflects the SU(3) symmetry. Highly de- 
generate "SU(3)-ferromagnetic" ground states at TJ% = 
are easily lifted either to give ferromagnetic states for 
U2 < or a spin-singlet for U2 > (when the boson 
number is even). 

The robustness of Theorems 1, 2, and 3 and their 
proofs are also impressive if one recalls that magnetic 
properties of the Fermi-Hubbard model depend crucially 
on lattice structures and electron numbers. Recall, for 
example, that Lieb's theorem [9|, [llj (which is a coun- 
terpart of Theorem 2) is valid only for the model on a 
bipartite lattice at half-filling (and its proof is quite in- 
genious). This is a manifestation of the fact that "ex- 
change interaction" in bosons are much more tractable 
than those in fermions. 

A remaining challenge is to take these advantages to 
give a concrete characterization of nontrivial magnetic 
structures for U2 > which are intrinsic to spinor bosons. 

It is a pleasure to thank Takuya Hirano, Tohru Koma, 
Fumihiko Nakano, Masahiro Takahashi, Akinori Tanaka, 
and Yuta Toga for valuable discussions. 
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